Conservation laws for voter-like models on directed networks 
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We study the voter model, under node and link update, and the related invasion process on a 
single strongly connected component of a directed network. We implement an analytical treatment 
in the thermodynamic limit using the heterogeneous mean field assumption. From the dynamical 
rules at the microscopic level, we find the equations for the evolution of the relative densities of 
nodes in a given state on heterogeneous networks with arbitrary degree distribution and degree- 
degree correlations. We prove that conserved quantities as weighted linear superpositions of spin 
states exist for all three processes and, for uncorrelated directed networks, we derive their specific 
expressions. We also discuss the time evolution of the relative densities that decay exponentially to 
a homogeneous stationary value given by the conserved quantity. The conservation laws obtained in 
the thermodynamic limit for a system that does not order in that limit determine the probabilities 
of reaching the absorbing state for a finite system. The contribution of each degree class to the 
conserved quantity is determined by a local property. Depending on the dynamics, the highest 
contribution is associated to influential nodes reaching a large number of outgoing neighbors, not 
too influenceable ones with a low number of incoming connections, or both at the same time. 



I. INTRODUCTION 

Conservation laws are intimately related to symmetries 
in the systems they hold for. They play an important 
role in the characterization and classification of differ- 
ent nonequilibrium processes of ordering dynamics. For 
example, in Kinetic Ising models one distinguishes be- 
tween Glauber (spin flip) and Kawasaki (spin exchange) 
dynamics. Kawasaki dynamics fulfils a microscopic con- 
servation law, such that the total magnetization is con- 
served in each individual dynamical step of a stochas- 
tic realization. This conservation law does not hold for 
Glauber As a consequence, the Glauber and Kawasaki 
dynamics give rise to different scaling laws for domain 
growth in coarsening processes pj, and they define dif- 
ferent nonequilibrium universality classes. 

In other types of nonequilibrium lattice models non- 
microscopic conservation laws are known to hold. They 
are statistical conservation laws in which the conserved 
quantity is an ensemble average defined over different re- 
alizations of the stochastic dynamics for the same distri- 
bution of initial conditions. Examples of such conserva- 
tion laws occur for the voter model 0, H[ or the invasion 
process 4]. In particular, the role of the conservation law 
of the magnetization and of the Zi symmetry (1 states) 
in the voter dynamics universality class has been stud- 
ied in detail in the critical dimension d = 2 of regular 
lattices [f|. The voter model is a paradigmatic model 
of consensus dynamics in the social context @, 0] or, in 
the biological context, of competition of plant species in 
ecological communities Q . In general, any Markov chain 
with at least two absorbing states reachable from all other 
configurations has a conserved quantity when averaged 
over the ensemble. Such a quantity determines the prob- 
ability to eventually reach a particular absorbing configu- 
ration in a finite system. In some cases, this conservation 
law is of rather trivial nature as in the zero temperature 



Ising Glauber dynamics where the magnetization sign is 
conserved. The voter model, the zero temperature Ising 
Glauber dynamics, and other related models of language 
evolution Q or population dynamics (Io| . belong to the 
class of models with two absorbing states whil e ep idemic 
spreading dynamics, like the contact process [JJJ or the 
Susceptible-Infected-Susceptible model [l2j , usually have 
a single absorbing state with no conservation law. 

While some of these questions have been studied for 
spin lattice models for a long time, conservation laws for 
dynamical processes on complex networks [l3|, [HI, HH, EH 
still remain a challenge. This issue has been considered 
for the voter model Q, Q or the invasion pro cess Q on 
undirected uncorrelated networks [l7], [Tl. 1 19L |20| . The 
link-update dynamics for the voter model has been found 
to conserve the global magnetization [2l| , while the node 
update dynamics [2lJ and the invasion process pre- 
serve a weighted global magnetization where the contri- 
bution of each spin is calibrated by some function of the 
degree of the corresponding node in the undirected net- 
work. Such ensemble average conservation laws charac- 
terize processes with two absorbing states accessible to 
the dynamics, that compete to maintain an active state 
in the thermodynamic limit. In finite networks, the con- 
served quantities give the probabilities of reaching the 
uniform states and so act as a bridge that enables some 
probabilistic predictive power of the final dynamical state 
based on information about the initial conditions. In ad- 
dition, different finite size dynamical scaling properties 
can be related to different conservation laws [2lj . 

Much less has been done exploring dynamical processes 
on directed networks, with the exception of the Ising 
model [221 ] and Boolean dynamics mainly applied to bi- 
ological problems [H|. However, interactions between 
pairs of elements are asymmetric in different systems in- 
cluding some social networks [24j, where social ties are 
perceived or implemented differently by the two individ- 
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uals forming the connected pairs. Directed network rep- 
resentations rather than undirected ones become more 
informative and adjusted to reality. In general, directed 
networks present characteristic large-scale connectivity 
structures, the so-called bow-tie architecture formed by 
a strongly connected component as a core structure and 
peripheral in- and out-components [25| . This organiza- 
tion, coupled to the initial condition of the dynamics 
running on top, have an impact both on the evolution 
of the processes and the final possible states of the sys- 
tems (2(| [53, [H|]. In the voter model, leaf nodes in the 
in-component never change their state thus sending an in- 
variable signal that can potentially propagate to the rest 
of the components of the system. This is closely related 
to phenomena such as the presence of zealots [2^, |3(| in 
undirected networks. Both input or output directional 
large-scale components and zealotry imply at the end an 
external forcing on the dynamical processes that prevents 
reaching one of the absorbing states even for a finite net- 
work. This is clearly illustrated by the evolution of dy- 
namical processes running on networks at the transition 
from a pure strongly connected component to a complete 
bow-tie structure. In an isolated and strongly connected 
component, the voter dynamics keeps an active dynam- 
ical state in the thermodynamic limit, but it leads to a 
consensus (absorbing state) in a finite network as it hap- 
pens on undirected networks. Thus, the appearance of 
an input component in the large-scale structure of the 
network prevents the system from reaching an absorbing 
state for random initial conditions [27l |. 

In this paper, we focus on dynamics of coupled two- 
state spin variables and consider conserved quantities 
that are weighted sums of the spin values. Specifically, we 
investigate the form of the conservation law for the voter 
model — under node and link update — and the invasion 
process in directed networks with arbitrary degree distri- 
bution and degree-degree correlations. The directionality 
of the interactions is therefore encoded in the topology. 
We restrict to a single strongly connected component so 
that the absorbing state can be reached in a finite sys- 
tem, what seems realistic for a number of densely con- 
nected real networks like the world trade web [31(. In 
Sec. II, we present a detailed study of the node update 
version of the voter model and implement an analytical 
treatment using the heterogeneous mean field assump- 
tion in the thermodynamic limit. From the dynamical 
rules at the microscopic level, we find the equations for 
the evolution of the relative densities of nodes in one of 
the two possible states on heterogeneous networks with 
arbitrary degree distribution and degree-degree correla- 
tions. In this case, we prove that a conserved quantity 
as a weighted linear superposition of spin states exists. 
In Sec. Ill, we discuss the node-update voter model in 
uncorrelated directed networks to derive analytical ex- 
pression for the conservation law and we also discuss the 
exponential decay of the relative densities to their homo- 
geneous stationary value, which is basically a function 
of the conserved quantity. We show how the conserved 



quantity determines the probability of reaching one of 
the two states in a finite network. In Sec. IV and V, 
we present the results of applying the same methodol- 
ogy to the voter model with link update and the invasion 
process, respectively. We conclude in Sec. VI with a sum- 
mary of results and open questions for future research. 



II. THE VOTER MODEL ON STRONGLY 
CONNECTED COMPONENTS 

In the voter model under node update (VM), each node 
of a network can exist in one of two possible states, 1 or 
46]. In a single dynamical event, a randomly selected 
node copies the state of one of its neighbors, also selected 
at random. The link update dynamics of the Voter model 
selects instead a link [21] . Time is increased by 1/N, so 
that the physical time is incremented by 1 after N of such 
events. On undirected networks, the node-update voter 
model conserves the ensemble average of a weighted mag- 
netization, where the contribution of each spin is multi- 
plied by the degree of the corresponding node. 

As defined above, the interactions in the voter dynam- 
ics are instantaneously asymmetric since the updates al- 
ways go in the same direction once the original node is 
chosen independently of the undirectionality of the sub- 
strate. Hence, the discussion of the voter model on di- 
rected networks comes out as a natural one, where the 
directionality of the interaction is decoupled from the dy- 
namics and encoded in the structure of the substrate. 
The straightforward generalization of the voter model on 
directed networks under node update consists of select- 
ing a node at random, and then assigning to it the state 
of one of its incoming neighbors, also chosen at random. 
We will discuss this dynamics next in this section and 
Sec. Ill, and the voter model with link update will be 
discussed later in Sec. IV. 



A. Directed networks 

The topological structure of directed networks is more 
complex than the one of undirected graphs. In purely 
directed networks, without bidirectional links, the edges 
are differentiated into incoming and outgoing, so that 
each vertex has two coexisting degrees k- lri and fc out , with 
total degree k = k in + k out . Hence, the degree distribu- 
tion for a directed network is a joint degree distribution 
P(k in , fc out ) = -P(k) of in- and out-degrees that in gen- 
eral may be correlated. We consider degree correlations 
Pi n (k'|k) and P out (k'|k), which respectively measure the 
probability to reach a vertex of degree k' leaving from a 
vertex of degree k using an incoming or outgoing edge of 
the source vertex, and are related through the following 
degree detailed balance condition [32] 

fc out P(k)P out (k'|k) = fc( n P(k')Pin(k|k'). (1) 
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This ensures that the network is closed and (fc; n ) = 
(fc out ). Apart from the prescribed degrees and two point 
correlations, networks are maximally random. 

At the macroscopic scale, the giant weakly connected 
component, i.e., the set of nodes that can communicate 
to each other when considering the links as undirected 
[II, [II, [II, [II, [I?]], becomes internally structured in three 
giant connected components, as well as other secondary 
structures such as tubes or tendrils, forming a bow-tie 
architecture [25j . The main component is the strongly 
connected component (SCC), a central core formed by 
the set of vertices that can be reached from each other 
following a directed path. The other two main compo- 
nents are peripheral components, the in component (IN) 
formed by all vertices from which the SCC is reachable by 
a directed path but that cannot be reached from there, 
and the out component (OUT) formed by all vertices 
that are reachable from the SCC by a directed path but 
cannot reach the SCC themselves. Percolation theory for 
purely directed networks was first developed for uncorre- 
cted networks [H, H3, [H, [H, 53] , and directed random 
networks with arbitrary two point degree correlations and 
bidirectional edges [32 1 . 

We restrict to networks forming a strongly connected 
component without peripheral components that would 
act on the SCC as sources of external forcing. We will 
see that within the strongly connected component, con- 
servation laws preserve weighted magnetizations, where 
the weights are dictated by the directed degrees. 



B. From microscopic dynamics to the drift 
equation under the heterogeneous mean held 
assumption 

To study the time evolution of the system, we consider 
the drift part of the Langevin equation for the density 
of nodes in one of the states from a microscopic descrip- 
tion of the evolution of single nodes' states applying a 
heterogeneous mean field approach. To our knowledge, 
this methodology, which allows us to deal with dynam- 
ical processes on complex networks, was first presented 
in Refs. Jill . |42| and recently used to study the contact 
process In [2(| EH , a homogeneous mean field pair 
approximation was instead developed. 

We focus on the microscopic state of nodes at some 
time t. Let s u (t),u — 1,...,N, be a stochastic bi- 
nary variable defined for each of the N nodes in the 
network which describes its state, or 1. The vector 
s(t) = {s u (t)},u = 1,...,N, completely defines the dy- 
namical state of the system at time t. Two more in- 
dependent binary stochastic variables [i(dt) and £ M are 
defined in order to model the transitions between states 
of single nodes in an iteration. After a time interval dt, 
the variable fi(dt) for a given node u takes the value 1 
or if u was chosen or not, respectively. In case node u 
was selected, then assumes the value 1 [0] if u copies a 
neighbor with state 1 [0] . We assume that the occurrence 



of events in the voter dynamics follows an independent 
Poisson process for each node, with constant rate A for all 
of them, which corresponds to a Montecarlo step. In the 
remainder we be set to A = 1 without loss of generality. 
Thus, /i(dt) and have probability distributions 

P(p(dt)) = dt5 ll{dt)A + (1 - dt)S Kdt)t0 , (2) 

P(£,u) = ^uAu.in^l + (1 - *«/fc«,in)^ u ,0, (3) 

where k u ^ n is the incoming degree of node u, and we 
have defined $ u (t) = a vu s v (t). The adjacency matrix 
{«„„} encodes the topological properties of the directed 
network. Element a vu has value one if there is a directed 
link from v to u and zero otherwise, so that <&«(£) stands 
for the number of state-one incoming neighbors of node u 
at time t. The matrix {a vu } is symmetric for undirected 
networks but for directed ones it is in general asymmetric. 

In terms of the above variables, the dynamical state 
s u (t) of node u after an increment of time dt is 



s u (t + dt) = ix{dt)£ u + (1 - n{dt))s u (t). 



(4) 



This equation, together with Eqs. ^ and (J3J, give the 
complete description of the evolution of the system, mak- 
ing the formalism general and applicable to any network 
structure. 

Although exact, this microscopic description is unman- 
ageable. In order to reduce the degrees of freedom, we 
apply a heterogeneous mean-field hypothesis [HI] so that 
nodes with the same degree k are assumed to be statisti- 
cally independent and equivalent and can be aggregated 
in the same degree class T(k) = T(ki n ,k out ). Properties 
are then defined for each degree class, that will be char- 
acterized by the relative density m^t), the ratio between 
the number of state-one nodes within class T(k) and its 
number of nodes N^, 



m k (i) 



;T(k) 



Su(t) 



(5) 



In the thermodynamic limit, the relative densities m k (i) 
can be considered as continuous variables. Their time 
evolution can be described by a Langevin equation [44j 
with drift and diffusion coefficients that are respectively 
given by the first and second infinitesimal moments of 
the stochastic variables m k (i). Those moments can be 
derived from the microscopic equation Eq. Q along with 
the definition in Eq. ([5]). In the thermodynamic limit, it 
is possible to prove that the diffusion term has a depen- 
dence \UfNk on the system size as for undirected net- 
works [20J, so that the drift term will dominate. It 
is given by the average value over all possible configura- 
tions of m-k(t + dt) conditioned to the state of the system 
at time i, 

(m k (t + dt)) mk(t) = m k (*) + A k (t)dt. (6) 
From the microscopic dynamics 

(s a (t + dt)) s{t) = s u (t) - dt [s u (t) - $ u (t)/k uM } , (7) 
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and summing this equation for all nodes in the degree 
class k and dividing by the number of nodes N^, we arrive 
to 



(m k (t + dt)) mk{t) = m k (t)-dt 



and from here to 
A k (t) = - 



K m ueT(k) 

(8) 



-Wk(i) 



1 1 



E 

ueT(k) 



(9) 



The adjacency matrix contained in can be coarse- 

grained as well, so that a differential equation for the 
relative densities can eventually be written. This coarse- 
graining restricts the validity of the equations to random 
complex networks (and not lattices), since we assume all 
nodes in the same degree class to be statistically inde- 
pendent. With these assumptions, 



weT(k) k' u6T(k') ueT(k) 



Q>vu,Sv (^) 



k' 



flk'kiVkiVk'mk'(t). 



(10) 



At this point, we restrict to directed networks organized 
at the large scale into a SCC without IN and OUT. This 
allows us to write 



Ok'k — 



^ ut P out (k|k') fc in P in (k'|k) 



N v 



(11) 



where Evk is the asymmetric matrix of the number of 
connections from the class of vertices of degree k' to the 
class of vertices of degree k, and where we have made use 
of the detailed balance condition Eq. (JTJ. 

Inserting these results into Eq. ©, we arrive to the 
equation for the evolution of the relative density in the 
degree class k of a purely directed correlated network 
(disregarding diffusion terms), 

^^ = -m k (i)+£Pin(k'|k)m k ,(i). (12) 

k' 

Let us recall that this result is valid for the ensemble 
of networks defined by the degree distribution P(k) and 
the degree correlations Pi„(k'|k) and P out (k'|k), but oth- 
erwise maximally random. Notice that big enough net- 
works present good statistical quality at the level of de- 
gree classes and are also well described by this equation. 
Finally, in the thermodynamic limit, the Langevin equa- 
tion loses its noise term because of the dependence on 
the system size and reduces to Eq. (TT2"1) . so that m,k(t) 
becomes a deterministic variable. Nevertheless, since the 
process is linear, Eq. ([12]) is always valid even for finite 
systems understanding that in this case the variables are 
averages over realizations of the process with the same 
distribution of initial conditions. 



C. Conserved quantity on directed networks with 
degree-degree correlations 

For correlated networks, mk(i) = J2 k , Pin(k'|k)m k '(i) 
in the stationary state and hence all relative densities are 
sfitangled through topological correlations. This equa- 
tion corresponds indeed to an eigenvector problem, since 
{m k (i)} can be thought as the eigenvector of the matrix 
{P; n (k'|k)} with eigenvalue one. 

We prove next that, within the heterogeneous mean 
field approach and for the correlated directed networks 
we are considering, there is a conserved quantity given 
as a linear superposition of the form u> — (p^m^t) . 
From Eq. (O, its evolution is given by 



du> 
~dt 



EI>kPn(k'|k)m k ,(i), (13) 



and imposing that dio/dt = 0, we obtain 

2J ^mk(t) = 2J £ -Pin(k|k')<,9 k /?7ik(i). 



k k' 



For each density 



^k = £Pi„(k|k> k ' 



(14) 



(15) 



This is an eigenvector equation that has a solution if the 
matrix {Pj n (k|k')} has an eigenvalue equal to one with 
{(/9k} the corresponding eigenvector. One can prove that 
this eigenvector with eigenvalue one exists by summing 
both sides of the previous equation over k. Using the nor- 
malization of the conditional probability ^ k Pj n (k|k') = 
1, one eventually arrives to a trivial identity [47:] . The 
fact that the coefficients ipu that modulate the contri- 
butions of the different ?7ik to the conserved weighted 
magnetization correspond to the entries of the eigenvec- 
tor of a certain characteristic matrix with eigenvalue one 
also applies to other similar dynamical processes, such as 
the link dynamics and the invasion process, as we will 
show. 

This proves that a conserved quantity of the form of a 
linear functional exists but, in general, it is not possible 
to derive its value without further specifying the form of 
the degree-degree correlations in the network. 



III. VOTER MODEL ON UNCORRELATED 

sees 

When two-point correlations are absent, the transition 
probabilities become independent of the degree of the 
source vertex. In this situation, 

P out (k'|k) = , P in (k'|k) = ^gp, (16) 

and using these expressions, Eq. ([12]) becomes 
Gfmk(f) 



dt 



-m k (i) + w ut! 



(17) 
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where we have defined 

^out = TT~T Y\ fcout^(k)m k (t). (18) 

\"W k 

Therefore, in the stationary state m k — uj out Vk and ui out 
is a conserved quantity in uncorrelated networks, which 
immediately follows from Eq. (|17[) . In general, it is not 
preserved in strongly connected components of directed 
networks with degree-degree correlations. This is in con- 
trast to undirected networks, where the conserved quan- 
tity li — {J2k kP(k)mk(t)) / (k) is preserved even in the 
correlated case and indeed for any structure [211 ]. Going 
back to the uncorrelated case, notice that the out-degree 
is the quantity that weights the contribution of the nodes 
to the conserved quantity. From a local perspective, what 
seems therefore important in the VM is to be able to in- 
fluence a large number of partners 

In uncorrelated networks, the convergence of the state- 
one relative densities to their stationary value can be 
easily computed. From Eq. (TTT)) . taking into account 
that w ou t is a conserved quantity and for a given initial 
condition mk(O), it is straightforward to arrive to the 
solution 

mk(f) = ^out + (mk(0) - Wout) e"*, (19) 

where we have substituted (fcj n ) by (/c ut)- Thus, all the 
densities decay exponentially fast to the stationary value 
TO k = w °ut an( i the relaxation time is for all of them 
equal and independent of the degrees. 

In the thermodynamic limit, the partially ordered sta- 
tionary state is stable, while finite-size fluctuations even- 
tually bring the system to one of the two possible una- 
nimity states. The probability P\ that the system ends 
up with all nodes in state one (m^ — 1, Vk) is given by 
the initial condition, that fixes the value of the conserved 
quantity at the beginning of the process. To see this, one 
takes into account that w ou t is an ensemble average con- 
served quantity of the form in Eq. (|18p . from which 

w«wt=-Pi- (20) 

This is in agreement with the fact that, in general, the 
Markov property of a stochastic process, if present, triv- 
ially ensures that the exit probability is a conserved quan- 
tity corresponding to a time-translation invariance. If 
the process has one absorbing state, the exit probability 
has a constant value one but, if the process has two or 
more absorbing barriers, the probability of reaching one 
of those is not trivial any more. 

It is also interesting to investigate what happens to 
the quantity Vi(t) = (J^ k h n P(k)mk(t)) /{km), which 
involves in-degree instead of out-degree. In the un- 
correlated case, and disregarding fluctuations, Vi(t) = 
(i»i(0) — w out )e _ * -I- w out , that is, in general Vi decays ex- 
ponentially fast to uj ut- The quantity Vi(0) depends on 
the initial condition. If this is homogeneous over degree 
classes, then Vi(0) = w Q ut and Vi{t) remains constant. 



In order to check the convergence of the sate one rel- 
ative densities to the conserved quantity, we have run 
numerical simulations of the voter model dynamics on a 
random uncorrelated network of size N = 10 5 , scale- free 
in-degree distribution with exponent 2.5 and exponential 
out-degree distribution. To obtain an initial state that is 
inhomogeneous in the densities m^, we have chosen an 
initial configuration in which half of the nodes with the 
lowest out-degree have state zero, and the other half have 
state one. In this way, initial densities nik(O) in classes 
with fc out lower than 4 were small or zero, while densities 
in classes with fc out larger than 4 were one. 

In Fig. [TJ we plot the average of the conserved quantity 
uj out and the densities for classes k = (fc; n , fe ut) = (2, 1), 
(4, 3) and (3, 9) vs time, over 100 independent realiza- 
tions starting from the same initial condition as men- 
tioned above. As predicted by the theory, we observe 
that (wout) stays constant over time, whereas the three 
densities converge to the average of the stationary value 
m k *, in a time of order 10. We note that, apart from finite 
size fluctuations, the convergence of the densities to m k 
happens for every realization. This can be seen in Fig. [21 
where we show the evolution of rn/2,1) and m(3,9) vs u> out 
in a single run. After a short transient, the densities and 
the conserved quantity start to evolve in a coupled man- 
ner (except from small deviations around the TOk = cj out 
line), they fluctuate from to 1 until they reach the ho- 
mogeneous zero-state. We also observe that fluctuations 
in TO( 3 g) are larger than in to^i), given that degree dis- 
tribution make the number of nodes in class (2,1) larger 
than in class (3, 9). 



IV. VOTER MODEL WITH LINK UPDATE 

The same assumptions and procedures apply to the 
link-update voter model and the invasion process. The 
link update (LU) dynamics selects first a directed con- 
nection, so that the node at the tail will always transmit 
its state to the neighbor at the head. 

The microscopic dynamics of the link-update voter 
model is described by 

s u {t + dt) = fi u (dt)£ u + (1 - n u (dt))s u (t), (21) 

where as for the voter dynamics £ u is given by Eq. Q 
and the binary variable fi u (dt) for the selection of a link 
has a probability distribution 

(dt),l + (1 — ^uan^)^„(dt),0- (22) 

A factor A/ (N (fc; n )) has been reabsorbed in the definition 
of dt. Proceeding as for the voter model (we skip the 
details), we arrive to the equation for the evolution of 
the relative densities for the different degree classes, 

^£ - -k in m k (t) + fc in ^P in (k'|k)TO k ,(t). (23) 
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FIG. 1: Time evolution of the conserved quantity o; ut (cir- 
cles) and the densities of state-one nodes in degree classes 
k = (fen, fcout) = (2, 1) (squares), (4,3) (diamonds) and (3,9) 
(triangles) , for the voter model dynamics. Curves correspond 
to averages over 100 realizations on a single random uncorre- 
cted network with TV = 10 nodes, scale-free in-degree distri- 
bution with exponent 2.5 and exponential out-degree distri- 
bution. While {lj out) remains roughly constant over time, the 
densities quickly decay to the stationary value (uwt)- The 
inset shows that the ratio between the densities and the con- 
served quantity is close to one during the entire evolution. 




out 



FIG. 2: Densities of state-one nodes m(a,i) and Wi(3,g) vs Wout 
in a single realization of the voter model dynamics on the same 
network of Fig. [1] The trajectories of classes (2, 1) and (3, 9) 
start at the positions (0.8,0) and (0.8, 1.0) respectively, then 
they quickly hit and move along the diagonal mk = w ut , until 
they reach the zero-state consensus point »?i(2,i) = m (3,9) = 0. 



Regarding the stationary state, the same result as for the 
voter model is found. The state-one relative densities be- 
have again as Wk(i) = J^ k , Pi n (k'|k)mk'(i), so that all 
the relative densities are entangled through topological 
correlations. We can once again prove, within the hetero- 
geneous mean field approach and for correlated strongly 
connected components, that a conserved quantity of the 



form m — ^ k (p^m^t) exists and is defined by the eigen- 
vector problem 



<^k = ^Pi n (k|k')<Pk', 



(24) 



where now (p-^ — fe n <pk- In general, it is not possible 
to derive these coefficients without further specifying the 
form of degree-degree correlations in the network. 
When two-point correlations are absent, 



dm\t(t) 
~dt 



-k in m k (t) + k in uj out (t). 



(25) 



In the stationary state, — uj out (t) Vk, but w out (i) is 
not a conserved quantity for the link update process as it 
was for the voter model. Instead, the conserved quantity 
is 

"'in / \ A, ir , 



£^P (k)mkW/( fc - 



(26) 



which follows from Eq. (|25)) . Compare this expression 
with that for the total magnetization in uncorrelated 
undirected networks w — ui = (J2k P(k)rrik(t)) / (k) 
which corresponds to the conserved quantity for those 
structures [21[. The dependence of the conserved 
weighted magnetization on the ratio between out- and 
in-degree for directed networks highlights the fact that in 
LU it is important to have both a high out-degree to be 
influential and at the same time to have a low in-degree 
not to be too influenceable. Notice that the ratio of the 
directed degrees is well defined since we are assuming that 
the network is organized at the macroscopic scale into a 
SCC without peripheral components all nodes having at 
least one incoming and one outgoing link. Finally, in fi- 
nite systems the probability of the state-one absorbing 
state is given by the conserved quantity, w ; = Pi, and 
so fixed by the initial condition. 

The derivation of how the state-one relative densities 
converge to their stationary value in uncorrelated net- 
works is more intricate than for the voter model, but we 
can make use of a quasi-stationary approximation [44] in 
order to solve Eq. (|25|) . exploiting the fact that w ; is the 
conserved quantity. In the stationary state to ou t = u i, 
and we approximate the equation by 



drrik(t) 
dt 



-k in m k {t) + fc in w„ 



(27) 



For a given initial condition mk(0), the solution is 

m k (i) = w oi (m k (0) - w oi )e- fci " t . (28) 

As in the voter model, all the densities decay exponen- 
tially fast to the stationary value w i, but in contrast 
not all the densities decay with the same velocity, which 
depends on the in-degree. Higher in-degree classes have 
smaller relaxation times and decay faster than lower ones, 
but the transient is always faster as compared to the VM. 
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TABLE I: Conserved quantities for voter-like models in strongly connected components of directed networks. 1st column, 
existence of conserved quantity for correlated networks; 2nd column, conserved quantity for uncorrelated networks; 3rd column, 
stationary values for the relative densities; 4th column, density decay. 



VM 
LU 

IP 



W„ 



m k 



m k 



Wont = p^-y E k fcout-P(k)m k (£) aw m k + (m k (0) - m k )e~ 

" oi = <fcout/fci„) ' 



3 Woi = i E k fcout/fc in -P(k)m k (t) w oi m k * + (m k (0) - m k *)e fei "* 



(l/fc in > ^ k fc in 



E k ^P(k)m k (t) 



Win m k ' + (m k (0) - m k *)e < fc in> 



V. INVASION PROCESS 

The invasion process (IP) picks nodes at random that 
export their state to a randomly chosen outgoing neigh- 
bor. A certain node u will update its state in a passive 
form only when one of its incoming neighbors v is se- 
lected as the first node in one iteration of the dynamics 
and then v chooses u among all its outgoing neighbors to 
transmit it its state. In this situation, it is more conve- 
nient to work with the probability of node u undergoing 
a state update with final state 1, £,u\ and the probabil- 
ity of node u undergoing a state update with final state 
0, £i . The probability distributions of these dichotomic 
stochastic variables are 

m { u 1] ) = KdtS^ tl + (1-^)6^^ (29) 

P(Z { u 0) ) = *ldt5 w 1 + (1 - $ldt)S m , (30) 



with 



^ ^ Q>vu s v (^) / k V) out j (31) 

V 

y,QW(l - s v (t))/k 

v,out 

(32) 



and the parameter A of the Poisson process for the hap- 
pening of events reabsorbed in dt. Using these expres- 
sions, the dynamics is described at the microscopic scale 
by 

s u (t + dt) = eHdt)(l-^(dt)) 

+ il-£\dt))(l-&°\dt))s u (t). (33) 

Following the same methodology as for the voter 
model, the drift equations for the relative densities in 
the different degree classes read 



dm k (i) 
dt 



= knY, p!-Pin(k'|k)(m k ,(i) - m k (t)). (34) 



The existence of a conserved quantity u> = E k ^k^kW 
in the correlated case is governed by the eigenvalue prob- 
lem 



. ^ P in (k|k')/fc out 



-f E k »^(k"|k')/fc. 



out 



(35) 



where <^ k = ^k^in Ek" • p m( k "l k )/ fc out- Summing both 
sides of this equation over k, one arrives once more to 
a trivial identity and so a conserved quantity exists in 
general on networks with degree-degree correlations. As 
we see next, we can be more specific on uncorrelated 
networks, for which Eq. ((34]) reduces to 



dm k (t) 
dt 



k lr 



-(m(t) - m k (t)), 



(36) 



where m[t) = Ek -P( k ) m k(0 is the total density of state- 
one nodes in the network. 

In the stationary state, m k (t) = m(t) Vk, but here 
m(t) is not a conserved quantity for the IP in uncorre- 
lated directed networks. Instead, the conserved quantity 
is 



>(*) 



m k (t) 



1 



&in / \ fe : . 

= Er p ( k KW/(r 



(37) 



In finite systems, the probability of the state-one absorb- 
ing state is given by this conserved quantity, oj- m = P\, 
and is therefore fixed by the initial condition. The de- 
pendence of the weights on the inverse of the in degree 
implies that those nodes with low in-degree, so less influ- 
enceable, have the highest contribution and control the 
process. This dependence on the in degree is analogous 
to the dependence on the degree of the conserved quan- 
tity w = uj = (Efc l/kP(k)mk(t)) /(k) in uncorrelated 
undirected networks [191 ]. 

After a transient, m(t) reaches the value uji n , so that 
the stationary values of the relative densities are m k (i) = 
Win Vk. This result tells us that all the densities become 
independent of k and reach the same stationary value, as 
in the previous processes. 

The derivation of how the state-one relative densities 
converge to their stationary value in uncorrelated net- 
works is more intricate than for the voter model, but like 
for the link update we can make use of a quasi-stationary 
approximation [44| in order to solve Eq. (|36p . Substitut- 
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ing into Eq. (|36p that in the stationary state m(t) = Lu- m , 

For a given initial condition mk(0), the solution is 

m k (i) = Win + (m k (0) - w in ) e . (39) 

All the densities decay exponentially to the stationary 
value u>- m . Higher in-degree classes decay faster than 
lower ones with a relaxation time that is proportional 
to the inverse of the in-degree, as is the case for LU. Due 
to the average degree in the relaxation time, however, 
transients are generally slower in the IP than in the LU. 
When compared with the VM, the IP dynamics exhibits 
a slower transient for degree classes with in-degree be- 
low average while those with in-degree above the average 
converge faster to the stationary state. 



VI. CONCLUSIONS 

We have introduced an analytical formalism from mi- 
croscopic dynamics to show that three different nonequi- 
librium dynamical models with two-absorbing states run- 
ning on strongly connected components of directed net- 
works with heterogeneous degrees and degree-degree cor- 
relations have associated ensemble average conservation 
laws. These conservation laws have been fully determined 
when degree-degree correlations are absent. The exis- 
tence of ensemble average conservation laws is a general 
characteristic of Markov processes with two or more ab- 
sorbing states. 

The constraints imposed on the dynamics by the con- 
servation laws lead to interesting and nontrivial behavior. 
From a practical point of view, they are related to the 
stationary values and the characteristic relaxation times 
of the relative densities of nodes in state one in each de- 
gree class and, in finite systems, gives the probabilities 
of reaching the two possible absorbing states. In this 
sense, the conservation laws obtained in the thermody- 
namic limit for a system that does not order in that limit 
(i.e. does not reach the absorbing state) determine the 
probabilities of reaching each absorbing state for a finite 
system. The contribution of each node to he conserved 
global weighted magnetization is always a specific func- 
tion of the directed degrees. In the case of the VM, the 
out-degree is the weight that controls the importance of 
the node as a measure of its influence, while in the IP 



it is the inverse of the in-degree, and in the LU it is the 
ratio between out and in-degree. In all cases, the con- 
served quantities are determined by local properties that 
encode the importance of each node in the network. De- 
pending on the dynamics, what seems important from 
a local perspective is to be influential reaching a large 
number of neighbors, or not to be too influenceable, with 
a low number of incoming connections, or both at the 
same time. 

From a broad perspective, these studies help in the un- 
derstanding of how the rich structure of real systems af- 
fects the dynamical processes that run on top. However, 
many questions still remain unsolved. In which specific 
way do degree correlations alter the results for uncorre- 
cted networks? How is the diffusive fluctuations regime 
in SCCs of finite directed networks? Is the finite size 
scaling of consensus times the same as in undirected net- 
works? On the other hand, it seems realistic to restrict to 
SCCs for a number of densely connected systems, like for 
instance the world trade web [3l[ , but in sparse directed 
networks the whole structure of core and peripheral com- 
ponents should be taken into account. Numerical simula- 
tions in some specific model networks [27j show that the 
appearance of an input component seems to prevent the 
system, even if finite, from reaching an absorbing state 
for specific initial conditions. How does the complete 
structure of a directed network couples to the initial con- 
ditions of the dynamics to induce the presence of zealots 
and how do they affect in quantitative terms the behavior 
of the whole system still needs further research. 

During the final completion of this work, we became 
aware of a recent preprint [45j discussing the fixation 
probabilities of mutants for Voter-like dynamics on di- 
rected networks. Since there exists a direct relation be- 
tween fixation probabilities of mutants and exit proba- 
bilities, and so conserved quantities, some of the results 
derived in that paper -without reference to conservation 
laws- concerning the dependence on the directed degrees 
are in correspondence to some of our results on uncorre- 
cted strongly connected components. 
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